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DYNAMICAL CORRELATION FUNCTIONS FOR PRODUCTS OF 

RANDOM MATRICES 


EUGENE STRAHOV 

Abstract. We introduce and study a family of random processes with a discrete time 
related to products of random matrices. Such processes are formed by singular values 
of random matrix products, and the number of factors in a random matrix product 
plays a role of a discrete time. We consider in detail the case when the (squared) sin¬ 
gular values of the initial random matrix form a polynomial ensemble, and the initial 
random matrix is multiplied by standard complex Gaussian matrices. In this case we 
show that the random process is a discrete-time determinantal point process. For three 
special cases (the case when the initial random matrix is a standard complex Gaussian 
matrix, the case when it is a truncated unitary matrix, or the case when it is a stan¬ 
dard complex Gaussian matrix with a source) we compute the dynamical correlations 
functions explicitly, and find the hard edge scaling limits of the correlation kernels. The 
proofs rely on the Eynard-Mehta theorem, and on contour integral representations for 
the correlation kernels suitable for an asymptotic analysis. 


1. Introduction 

One of the most recent achievements in Random Matrix Theory is the discovery of the 
exact solvability of ensembles dehned by prodncts of random independent matrices. In 
particular, it was shown in Akemann and Burda [Ij, Akemann, Ipsen, and Kieburg [5] 
that both eigenvalues and singular values of products of complex Ginibre matrices form 
determinantal point processes, and that the kernels of these processes can be computed 
explicitly in terms of special functions. Different scaling limits of these determinantal 
point processes and of their extensions were investigated in the papers by Kuijlaars and 
Zhang [21], Kuijlaars and Stivigny [23], Liu, Wang, and Zhang [25], Forrester and Liu 
ini- We refer the reader to the paper by Akemann and Ipsen ]6] for a survey concerning 
exactly solvable ensembles dehned by products of random matrices. 

Such ensembles arise in the theory of wireless communication, in particular in the 
problem about transmission of vector-valued signals through a sequence of clusters of 
scatterers, see Muller [2B]- In this problem the transmission of a vector-valued signal 
between two adjacent clusters of scatterers is modeled by a random matrix, and the 
transmission of a vector-valued signal through a sequence of clusters of scatterers is 
modeled by the corresponding product of random independent matrices. Exact solvabil¬ 
ity of the relevant random matrix ensemble enables to give formulae for the expectation 
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of the mutual information, which is an important quantity in wireless telecommunica¬ 
tion, see Akemann, Ipsen, and Kieburg [S], and Wei, Zheng, Tirkkonen, and Hamalainen 

ra- 

Considering the model described above we note that the ordinal index of a cluster in 
the sequence can be associated with a discrete time, and the number of factors in the 
relevant random matrix product can be identihed with a discrete time as well. Motivated 
by this observation, the present paper deals with random processes with a discrete time 
constructed in terms of products of random matrices as follows. Recall that complex 
Ginibre matrices are those whose entries are i.i.d standard complex Gaussian variables. 
Assume that at time A; = 1 we are given a random matrix Gi of size Ni x Nq, where 
> Nq. Assume further that squared singular values of Gi form a polynomial ensemble, 
see Section [2A] for a precise dehnition. Set W(l) = Gi. At time k = 2 the matrix Gi is 
multiplied (from the left) by a new matrix G 2 , where G 2 is a complex Ginibre matrix of 
size N 2 XN 1 (where N 2 > Nq) independent on Gi. At time k = 3 the matrix X(2) = G 2 G 1 
is again multiplied from the left by a new complex Ginibre matrix G 3 of size N 3 x N 2 
independent on Gi, G 2 . Gontinuing in this way we obtain a sequence of product matrices 
{A(/c)}^^ dehned by X{k) = Gk ■ ■ ■ G 2 G 1 , where k = l,...,m, the squared singular 
values of Gi form a polynomial ensemble, and G 2 ,.. ■,Gm are independent matrices with 
i.i.d standard complex Gaussian entries. 

Note that for each k = 1,... ,m the matrices X*{k)X{k) are of the same size Nq x Nq. 
We wish to describe the evolution of the eigenvalues of X*{k)X{k). In other words, the 
problem is to describe the evolution of the (squared) singular values of X{k). Glearly, 
this evolution can be thought in terms of non-intersecting polygon paths parameterized 
by the (discrete) time parameter k. 

Denote by the jth largest eigenvalue of X{k)*X{k). If matrices Gk are independent, 
G 2 , ..., Gm are complex Ginibre matrices, and squared singular values of Gi form a 
polynomial ensemble, the conhguration (A;, x^) is shown to be a determinantal point 
process on {1,... ,m} x M>o. This is the main result of the present paper. For three 
special cases (when Gi is a complex Ginibre matrix, when Gi is a truncated unitary 
matrix, and when Gi is a complex Ginibre matrix with a source) we compute correlation 
kernels of the dynamical correlation functions explicitly, and hnd their hard edge scaling 
limits. The proofs rely on the Eynard-Mehta theorem [15], and on contour integral 
representations of the correlation kernels suitable for an asymptotic analysis. 

Let us mention that recent works on products of random matrices (in particular, the 
works mentioned above) deal with the situation when the number of matrices in products 
under consideration is hxed. In this article we propose to consider sequences of products 
of random matrices instead of hxed products. The main idea behind this work is rather 
similar to that of Johansson and Nordenstam [2T|. Johansson and Nordenstam [2T| start 
with an inhnite random matrix picked from the Gaussian Unitary Ensemble, and consider 
the eigenvalues of its minors. They show that all such eigenvalues form a determinantal 
process on N x R. Such determinantal processes are called the minor processes. Here 
we can start from a matrix product Gm ■ ■ - Gi , and consider squared singular values of 
sub-products ... Gi, A; = 1,..., m. As a result we arrive to a family of determinantal 
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processes on {1,... ,m} x M>o, which can be understood as time-dependent extensions 
of the processes studied previously. 

The minor processes and their generalizations are of great interest in Random Ma¬ 
trix Theory, and are closely related to random tiling models, and to certain percolation 
models on the Z"*" x lattice. In particular, it was observed by Johansson and Nor- 
denstam [21] (see also Okounkov and Reshetikhin [29]) that the spectra of the principal 
minors of a GUE-matrix behave as domino tilings of large size Aztec diamonds. We 
refer the reader to Adler, van Moerbeke, Wang [I], Adler, Chhita, Johansson, and van 
Moerbeke [2], Adler and van Moerbeke [5] and to references therein for a description 
of these connections, and for recent results in this area of research. Many others time- 
dependent deteriminantal point processes appeared over the past decade in the context 
of growth models, and were applied to different problems in numerative combinatorics 
and statistical physics, see, for example, Borodin and Gorin [10], Borodin and Ferrari 
[ 11 ] 

The novelty of the present paper is in the very observation that products of random 
matrices also lead to time-dependent determinantal point processes. These processes 
are rather different from those formed by eigenvalues of minors, have different scaling 
limits, and were not considered in the literature, to the best of author’s knowledge. One 
should expect that time-dependent determinantal processes originated from products of 
random matrices lead to interesting models of statistical mechanics (similar to the minor 
processes), but we leave this question for a future research. On the technical level the 
main contribution of the present paper is in the demonstration of a relation with the 
Eynard-Mehta theorem, and of its usefulness for derivations of correlation functions. 
This relation not only give new results (explicit formulae for the dynamical correlation 
functions), but also provides a new and a unihed point of view on the point processes 
formed by the squared singular values of random matrices. 

The paper is organized as follows. In Section |2] we state the main results obtained 
in this work. In particular. Theorem 12.11 of Section 12.11 gives a general formula for 
the correlation kernel of the determinant process under consideration. Sections 12.1112.41 
provide contour integral representations and hard edge scaling limits of the correlation 
kernels corresponding to three special initial conditions mentioned above. The rest of 
the paper is devoted to proofs of these results. 

Acknowledgements. I am deeply grateful to Dong Wang for a number of remarks and 
suggestions, and for bringing the preprint [13] to my attention. 


2. Main results 

2.1. Ginibre product process with polynomial ensemble initial conditions. Let 

Gi,.. .,Gm be random matrices, and assume that each matrix Gk, fc G {1,... ,m}, is of 
size Nk X Nk-i. Set 


X{k) = Gk ■ ■ ■ Gi, kE{l,...,m}. 
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If n = Nf), then for each fc, fc G {1,..., m}, X{k)*X{k) are random matrices of the same 
size n X n. Furthermore, set 


- ^0, j e , 

and assume that these parameters are positive. Denote by x’j the jth largest eigenvalue 
of X{k)*X{k). The conhguration of all these eigenvalues, 

{{k,x’'j)\k = 1,.. .,m; j = 1,.. .,n} , 

form a point process on {1,..., m} xM>o- The correlation functions of this point processjil 

gi : ({1,... ,m} X R>o)* -)■ M, 

will be called the dynamical correlation functions for the products of random matrices 
Gi,.. .,Gm- Also, we say that the random matrix Gi determines the initial conditions of 
the product matrix process under considerations. The random matrices G 2 , • • •, Gm will 
be called the transition matrices. 

In what follows we always assume that the matrices Gi,, Gm are independent, and 
the transition matrices G 2 ,..., Gm are complex Ginihre matrices. We consider the case 
when the squared singular values (xj,..., xlf) of the matrix Gi form a polynomial ensem¬ 
ble in the sense of Kuijlaars and Stivigny [2Sj. This means that the vector (xj,... ,x^) 
has densitjH 

n 

(2.2) const -A (xj,..., xl) det [/i(xj)]”^.^^ JJ 

where A (xj,... ,x\) = H (xl — x]) is the Vandermonde determinant. We assume 

that given functions /i, ..., /„ are such that expression fl2.2p can be understood as a 
probability density. We obtain a point process on {1,..., mj x R>o. Let us refer to this 
process as to the Ginihre product process with polynomial ensemble initial conditions. 
The density of this product process is that of the vector 

£= P‘,..., i”) 

where x^ = (x},... ,xl) are the squared singular values of Gi, and x^ = (x^,... ,x))), 
r G {2,...,m}, are the squared singular values of Gr ... Gi. 


^We refer the reader to Johansson for a general discussion of correlation functions for point 
processes. 

^Kuijlaars and Stivigny [23] defined the polynomial ensemble by the expression 
(2.1) const-A (xi,.. .,Xn)det . 


For a computational convenience, we use equation (ESI) instead of equation (EH). 
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To present our result for the Ginibre product process with polynomial ensemble initial 
conditions introduce the Meijer G-function (see e.g. Luke [2S]) 

(2.3) 

A 1 f nr=ir(&.-5)n;=ir(i-a,+5) 

) 27r^ J r(l - 6, + .) r(a, - ^ 

c 

Here 0<m<q, 0<n<p, and an empty product is interpreted as unity. The contour 
of integration C depends on the location of the poles of the Gamma functions, and we 
refer to the NIST handbook 1271 for details on the different possibilities. 


^m,n I 

I bi, 


0 ^2, 

h, 


Theorem 2.1. The Ginibre product process with polynomial ensemble initial conditions 
is a determinantal point process on {1,..., m} xR>o- Its correlation kernel can be written 
as 


Kn,m{r, x\ s,y)^- 


(2.4) 


X 


^r+1? • • • ? 


^s>r 


ij=l 


In the formula just written above the function (j)r,m+i{x,i) is defined as follows. Ifm = 1, 
then (f)r,m+i{.xG) = Form >2 we set 

. _ J +i).. .r(z/m+ z), r = l,...,m-l, 

(Pr,m+l{X,l) - ^ r = m. 

The matrix A = in the formula for the correlation kernel is given by 


(2.5) aij = T{iy 2 +j).. .T{iy^+j) / fi{t)e ^dt, i,je{l,...,n}, 


and the functions 0o,s o,re given by 

. fMv^^e-y, 

(2.6) (j)Q^^{j^y) = 


5 = 1, 






^ ^ dt, s = 2,... ,m. 


If the squared singular values of Gi form a polynomial ensemble, and G 2 ,..., Gm are 
independent Ginibre matrices, then Theorem 12.11 implies that the dynamical correlation 
functions pi for the products of random matrices Gi,.. .,Gm have the determinantal form, 

(2.7) gi (si, xi;...;si, xi) = det {Kn,m{si, xp, sj, , / = 1, 2,.... 

We note that Kn,m{p, x; s, y) has a structure similar to that of the correlation kernel for 
a continuous version of the Schur process, see Borodin and Peche [T^ Thm. 3]. Both 
Theorem 12.11 and Borodin and Peche [121 Thm. 3] can be proved using the Eynard-Mehta 
Theorem |15j . 
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Once the formula for the correlation kernel is derived many probabilistic quantities of 
interest can be be computed. In the process interpretation, the evolution of the squared 
singular values oi X{k) = ■ ... ■ Gi can be thought of as consisting of n paths pa¬ 

rameterized by a discrete time. Given 1 < /i < ... < /^ < m, and subsets Oi, ..., 
Vtk of M>o, one quantity of interest is the probability that for all i, i & {1,..., k}, no 
path passes through Oj at time /*. It follows from Theorem 12.11 that this probability 
is expressible as a Fredholm determinant dehned by the kernel Kn^mir, x; s,y). Such 
Fredholm determinants can be studied by different methods, in particular in terms of 
differential equations, see Tracy and Widom [21] for partial differential equations derived 
in the context of Dyson processes, and Strahov [20] for systems of ordinary differential 
equations derived for singular values of products of independent Ginibre matrices. Also, 
paper by Bornemann |8] gives a systematic numerical treatment of Fredholm determi¬ 
nants arising in the context of Random Matrix Theory, and the methods of this paper 
might be applicable for the Fredholm determinant dehned by the kernel x] s, y). 

Theorem 12.11 gives a general expression for the correlation kernel, which works for any 
sequence of functions /i, ..., /„, for which fl2.2p is a probability density. For some special 
cases we can provide more explicit formulae for the correlation kernels, in particular we 
can derive contour integral representations suitable for a subsequent asymptotic analysis. 


2.2. The initial matrix is a complex Ginibre matrix. Assume that Gi is a complex 
Ginibre matrix. Then the density of the squared singular values of Gi is given by equation 
fl2.2p . where the functions /i, ..., /„ are dehned by 

(2.8) fi{x) = i G {1,... ,n} . 


Theorem 2.2. When the initial matrix Gi is a complex Ginibre matrix, the Ginibre 
product process is a determinantal point process on {l,...,m} x M>o. Its correlation 
kernel, Kn^rn{'f'iX',s,y), can be written ajl 


(2.9) A'n,m(G X] s, y) = —Gg ,0° 


X 


^r+1? • • • ? 


n—1 


^ J ls>r “1“ ^ 

p=0 


where the functions Pj,,p o^nd Qs,p are defined by 


( 2 . 10 ) 


Pr,pix') 


T{t-p) 


27ri/ nLor(t + i^i + l) 


x^dt, 


and 

( 2 . 11 ) 


nc+ioo 


Qs,p{x) 


ULo'^iu + ufi _ 


y '^du, c > 0. 


T{u-p) 

Here Sp is a closed contour that encircles 0,1,... ,p once in the positive direction, and 
such that Ret > —1 for t G Sp. 


^Note that t'o = 0. This fact will be used throughout the paper. 
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Theorem 12.21 enables us to derive a double integral formula for the correlation kernel. 
Namely, the following Proposition holds true 


Proposition 2.3. We have 

Kn,m{r,x-s,y) = ( 

X ’ \ T'r+l 

1 r , n;=or(w + ^^, + l)r(t-n + l)a:V“"' 

+ (2777)2 J ll^ Tit + uj + l)T{u-n + l) t-u ^ 

1 • V 

— ^—lOO ^ 

where T, is a closed contour going around 0,1,... ,n in the positive direction and such 
that Ret > —^ for t G S. 

If r = s, then the determinant process defined by s,y) is that formed by 

squared singular values of products of r independent complex Ginibre matrices. In this 
case the first term in the right-hand side of equation fl2.12l) becomes identically equal 
to zero, and the second term coincides with the correlation kernel obtained by Kuijlaars 
and Zhang m, see Kuijlaars and Zhang [2H Prop. 5.1]. By the same argument as in 
Kuijlaars and Zhang [23] we find the scaling limit of „j(l, x; s, y) near the origin (hard 
edge). 



Proposition 2.4. Assume thatx, y take values in a compact subset o/M>o, and assume 
that the parameters ui, ..are fixed. We have 


(2.13) 

lim „ (r, s, -) 

n^oon \ n n/ 


^s—r,0 
^0,s-r 
X 


l^r+lj 


^s>r 


{ 2 Triy 




du ® dt 


n;^or(^ + Uj + 1) sin TTM x^y 
rifco + 1) sin vrt t — u 


— 77—200 


where r,s ^ {1,... ,m}, and Yi is a contour starting from -l-cxo in the upper half plane 
and returning to -|-cxd in the lower half plane, leaving — | on the left, and encircling 
{ 0 , 1 , 2 ,...}. 


Proposition 12.41 implies that the scaling limits of the dynamical correlation functions 
do exist at hard edge. These scaling limits characterize the asymptotic behavior of the 
point process under considerations as n —)■ oo. 


2.3. The initial matrix is a truncation of a random unitary matrix. Let U be 

an / X / Haar distributed unitary matrix, and let Gi be the {n + ui) x n upper left block 
of U. We assume that z/i > 0, and that I > 2n + ui. 
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Proposition 2.5. The density of the squared singular values (xj,... ,x^) of Gi can be 
written as 

n 

const - A (a:J,..., xy ^ {x]Y" (l - x])X(o,i)^)) 

i=l 


where 


Xio,i){x) 


1 , xe(0,l), 

0, otherwise. 


Proof See Jiang [191 Prop. 2.1], □ 

Thns the sqnared singnlar valnes (x^,..., x^ of Gi form a polynomial ensemble, and 
the fnnctions /i,..., in eqnation fl2.2p are given by 

(2.14) /i(x) = e"=x*"yi-x)'"2"“^iX(o,i)(x), ie n}. 


Theorem 2.6. When the initial matrix Gi is a truncation of a random unitary matrix (as 
defined above), the correlation kernel of the Ginibre product process is given by equation 
where the functions PrY^)j Qs,p{y) are defined by 


Pr^p{pP) 



Vit-p) 


Vjt + l - 2n + l+p) , 

■ ■ ■y* __ y \ «.T hZ/ 

nj=or(^+A +1) 


and 


Qs,p{y) 


c+ioo 

[ _ nj=or(» +a) _ 

2 TTi J T{u — p)T{u — 2 n + I + 1 + p) 


C —200 

Here r, s G {1,..., m}, c > 0, and the contour Sp is chosen in the same way as in the 
statement of Theorem \2.2\ 


As in the previons case, a donble contonr integral formnla for the correlation fnnction 
can be derived. 


Proposition 2.7. When the initial matrix Gi is a truncation of a random unitary ma¬ 
trix, the correlation kernel, Kn,m{r,x;s,y), can be written as 

(2.15) 

A',.,m(r, i; s, y) = — 

X 

— ^+ioo 

1 f r ^ n'=or(w + A + 1) r(t + 1 - n)r(t + l-n + l) x*y-^ 
(27ri)2 J nLor(^ + A + 1) + 1 -'R)r(M + Z - n + 1) u-t 

1 . y' -J 

— ■^—200 ^ 

where r, s G {1,..., m}, and S is defined as in Pro'position \2.3[ 

Formnla fl2.15p enables ns to hnd the hard edge scaling limit of the correlation kernel. 
After snitable scaling, we obtain the same limiting kernel as in the case where the initial 
matrix is a complex Ginibre matrix. 
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Proposition 2.8. Assume that the parameters vi, ..Vm are fixed, and that I grows at 
least as 2n. Then for x, y taking values in a compact subset o/M>o 


(2.16) 

lim 


(I 


-K„ 


n]n 


X 


; s. 


’ (/ — n)n’ ’ {I — n)n 


_ r,0 


^ 1 1 

I -*-s>r 

X 


iflni) 


— TT—"iOO 


, / nj=o r(M + L-J + 1) sin Tiu x^y 

du ® —z—-—^-, 

/ n,=o+ 1) t-u 

s 


where r, s G {1,..., m}, and S is chosen in the same way as in Proposition^^ 


We note that for r = s the formulae fl2.15l) and 02.161) reduces to those found by 
Kuijlaars and Stivigny 


2.4. The initial matrix is a Ginibre matrix with a source. Assume that the 
random matrix Gi determining the initial conditions of the Ginibre product process is 
a sum of a complex Ginibre matrix G and a nonrandom matrix Q. In this case we use 
terminology of Forrester and Liu mi and say that the initial conditions of the Ginibre 
product process are given by a complex Ginibre matrix with a source. 

Proposition 2.9. Assume that the initial matrix Gi is given by 

Gi = G + Q, 

where G is a complex Ginibre matrix of size {n + ui) x n, and Q is a nonrandom matrix 
of the same size whose sguared singular values are (gi, ..., g^)- Then the density of the 
sguared singular values (t];, ... ,xlf) of Gi is proportional to 

n 

, xl) det (oFi (z/i + 1; qiX^j))^ .^^ JJ e“GL 

j=i 

Proof. See Desrosiers and Forrester m Prop. 5], and Forrester [16], §11.6. □ 

Proposition 12.91 implies that the squared singular values of the initial matrix Gi = 
G + Q form a polynomial ensemble defined by equation (I2.2|l . where 

fi{x) = ’ i e {1,..., n} . 

According to Theorem 12. li the Ginibre product process with the initial matrix Gi = G+Q 
is again a determinantal point process. To give an explicit formula for the relevant 
correlation kernel we follow to Forrester and Liu na, and introduce new functions, 


A {xl,... 
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and d)s(g;|/). Namely, we set 

(2.17) 


{2myV{vi + 1 ) 


- / dwi... dwr Wi 
/=! 

71 7r 


.Wl 


X 1^1 + 1; 


XU 


Wl . . .Wr 


In the formnla above 71 , ..., 7 ^ are paths starting and ending at — cxd and encircling the 
origin once in the positive direction. In addition, set 


C+200 


®s(g;i/) = — / dzy ^id{q]z)r{u 2 +z).. .r{us +z), 


where c > —min (z/i, ..., s G {1,..., m}, and 


id{q;z) = 


* 0-^1 {^1 + 1 ; —Qt) 


Tiui + 1 ) 


Theorem 2.10. When the initial matrix is a a Ginibre matrix with a source, the Ginibre 
product process is a determinantal point process on {1,... ,m} x M>o. Its correlation 
kernel, Kn^myW]s,y), can be written as 


X 


Iln^miTi Xi ^iV') ^0,s—r 


^r+l 5 • • • 5 


^ 11 

I -•-s>r 
X 


(2.18) 


H- du dvu^^e ^~^'“'^r{u]x)^s{v]y) - TT 

2'xi II -?/ _ 91 J. J. 


u + qi 


u — v-^-^v + qi 

i=i 


0 c 

where C is a counterclockwise contour encircling —qi, ..., —q^ but not u. 


Theorem 12.101 is a time-dependent extension of Forrester and Lin [T^ Prop 1.1]. In¬ 
deed, once r = s Theorem 12.101 gives a correlation kernel for the determinantal process 
formed by sqnared singular values of the random matrix product Gr .. .G 2 {G + Q). We 
note that once Theorem 12.11 is established, the result of Theorem 12.101 can be derived 
using the same arguments as in Forrester and Liu im. 

Theorem 12.101 is a starting point for an asymptotic analysis. Because of similarity of 
formulae in Theorem 12.101 and in Forrester and Liu m Prop 1.1] such an asymptotic 
analysis is a repetition of arguments from Ref. HI. In particular, time-dependent 
extensions of the critical kernel, and of the deformed critical kernel of Forrester and Liu 
m Thms 1.2, 1.3] are given by Proposition 12.Ill below. 
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Proposition 2.11. (A. Subcritical regime) Set qi = ... = qn = bn, and assume that 
0 < 6 < 1. For X, y in a compact subset of (0, cxd), and for fixed parameters vi, .. Um 


lim .. Kn,m ( r, ,, ; s ^ 


n-^-oo — b)n 
1 




y- 1 1 

I J- 5 >r 
X 


— l+ioo 


(27rz) 


, I ,Jlj=o^iu + ^j+ ^)smnux^y ^ ^ 

du (p dt^=pp: -—-r- 

nj=0 + L'j + 1) sin Tit t — u 


_ ^5 —TjO 

X \ Vr+l, ■■■.Vs 


^s>r 


1 


^ 1,0 

0 ,r- + 2 \ 


UX G 


5+1,0 

0,5+2 


Uq,Ui,...,Us 


uy ) du, 


where r,s E {1,..., m}, and S is chosen as in Proyosition \2.A 
(B. Critical regime) Set 

qi = ... = qn = n{l- rj^Y^. 

Then for x, y in a compact subset of (0, cx)), and for r in a compact subset o/M we have 


1 - ^ ir ( ^ . y \ 1 r ,0 

lim -^Kji Yn I ^5 ^7 s, 1 Gq 
n^oo^/n \ \/n \/n / X ’ \ Z^r+l, • • • , I^s 

19 19 

If f /M\ 1^1 e"™“2“+^^+2^^ 

- du dv {— ] - 

27 rzy y* \vJ u-v 

0 


X 


'-5>r 


^ '-^0,r+2 


UX 1 ^ 0+^2 


t'O, 


7 • • • 7 ^5 


uy 


^0 7 Z^l, . . . 7 

(C. Deformed critical regime) Set 

qj = \/naj, j = 1 ,... ,1, and qk = n(l — rjk = I + 1,... ,n, 

and assume that ai, ..., ai are strictly positive. Then for x, y in a compact subset of 
( 0 , CX)), and for t in a compact subset o/M we have 


]• ^ T/' f X y \ ^'■)' 

hm —=Kn,m r, —; s, ^ = -Gq 1 

n^oo Jn \ \/n \ n j x 


5—r,0 
r 


19 19 

- du dv { — ) - 

27rzi i* \v/ u-v 

0 


^r+l7 • • • 7 


y 11 

I J- 5 >r 
X 


i 


X 11 


n y+fTj 

?=i y 


- Z/q , , . . . , Z/ 7 - 


UX I Gntio 


0,5+2 


Z/Q, Ui,...,Us 


uy 
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Finally, let us mention few directions for a further research. Theorem 2.1 in a very 
recent paper by Clayes, Kuijlaars and Wang [13] suggests that similar methods (as in 
this paper) can be applied to processes formed by eigenvalues of certain sums of random 
matrices. Also, paper by Clayes, Kuijlaars and Wang [13] gives a general formula for the 
correlation kernel of the squared singular values of the product Gr ■ ■ ■ G 2 C 1 , where Gi 
can be any random matrices whose squared singular values form a polynomial ensemble, 
while G 2 , ■ ■ ■, Gr are independent complex Ginibre matrices, see Corollary 2.16 in Clayes, 
Kuijlaars and Wang [T3]. A question arises whether there is a time-dependent extension 
of this formula. In addition, the methods of the present paper can be extended to matrix 
products formed by truncated unitary matrices, and time-dependent extensions of results 
in Kieburg, Kuijlaars, and Stivigny [22], Forrester and Liu im can be obtained. 

The idea to investigate the dynamical correlation functions comes naturally from con¬ 
sidering the signal model of the MIMO system described in Wei, Zheng, Tirkkonen, and 
Hamalainen [32] • In this model the products of independent complex Ginibre matrices 
arise, and the results of Section 12.21 seem to be relevant in the MIMO communication 
setting. The author is not aware whether other models considered in this paper (i.e. the 
products of truncated random matrices, and the products of random matrices with a 
source) are useful for the theory of wireless communication. 

3. Measures given by products of determinants and the Eynard-Mehta 


THEOREM 


The aim of this Section is to recall one of mechanisms which gives rise to determinantal 
point processes, namely the case when the measure dehning a point process is given by 
products of determinants. We follow Johansson [2D], Section 2.3, and use the same 
notation and terminology as in Johansson’s paper (see also Eynard and Mehta [T3] . 
Borodin [S], Tracy and Widom [31]). 

Let n, m > 1 be two hxed natural numbers, and let Xq, Xm+i be two given sets. Let 
X be a complete separable metric space, and consider a probability measure on (X”)™ 
given by 



(3.1) 


m—1 



r=l 


In the formula just written above is fhe normalization constant, the functions 

ipr^r+i :XxX^C, r = l,...,m — 1 are given transition functions, 0o,i • ^0 X X —)■ C is 
a given initial transition function, and 4>m,m+i ■ X x X^+i —?• C is a given hnal transition 
function. Also, 



the vectors 



.0 

1 ) ■ 



-n ^m+1 
0 5 ^ 



.m+1 


n 


m+1’) 


n 
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are fixed initial and final vectors, and 

m n 

r=lj=l 

Here /i is a given Borel measure on X. Given two transition functions (j) and "0 set 


J X 

The following statement is often referred as the Eynard-Mehta theorem. 


Proposition 3.1. The probability measure pn,m{x)dfi{x) defines a determinantal point 
process on ,m} x X. The correlation kernel of this determinantal point process, 

x; s, y) (where r, s G {1, ..., m}, and x,y E X), is given by the formula 

n 

(3.2) Kn,m{r, x; S, y) = y) + i,j V)' 


The functions (j)r,s, and the matrix A = (ajj) (where i,j = 1 ,... ,n) are defined in terms 
of transition functions as follows 


(3.3) 


fir,six, y) 


((br¬ 


and 


(3.4) Oij = fio^m+i{x°,xf^^). 

Proof. See Johansson [20], Section 2.3. □ 


4. Proof of Theorem 12.11 


Our starting point is Proposition 14.11 which gives the distribution of squared singular 
values of a complex Ginibre matrix multiplied by a constant matrix from the right. 


Proposition 4.1. Let G be a complex Ginibre random matrix of size (n + z/) x I, and 
let Q be a non-random matrix of size I x n. Let qi,... ,qn be the squared singular values 
of Q, and assume that qj is not equal to zero for each j, j = 1,... ,n. IfY = GQ, then 
the squared singular values yi, ..., yn of Y have a joint probability density function on 
[0, cx))” proportional to 


A(i/i,...,i/n) , , 

UTT-L clet 

A(gi,...,g„) 


Vi 


e 


lu 

y 


n 


Li=i 


Proof. See Kuijlaars and Stivigny [23], Section 2. 


□ 


Recall that the Ginibre product process with polynomial ensemble initial conditions 
is formed by singular values of the products Gk ■ ■ ■ G 2 G 1 , k G {1,..., m}. It is assumed 
that the density of the squared singular values of Gi is given by equation fl2.2p . and that 
Gk, • • G 2 are independent complex Ginibre matrices. Proposition 14.21 shows that the 
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density of this point process can be written as a prodnct of determinants as in eqnation 

Proposition 4.2. Set Xq = {1, ... ,n}, X^+i = {!) • • • X = M>o. Let fi be the 

Lebesgue measure on M>o. Assume that the initial transition function, 

00,1 • , n} X R>o 

is defined by 

fioAfix) = fiix), 

and that the final transition function, 

0m,m+l • I^>0 X {1, . . . , nj- ^ C, 

is defined by 

fim,m+li,X, fi X 

In addition, assume that the transition functions 4>r,r+i '■ ®>o x M>o —j-C, r = — 

1, are defined by 

r = 1, 

(j)r,r+iix,y) = { , r = 2, 


X 


Pm — 1 — Pm — 1 /> f Pm ^ 


y’^rne a, ^ T = TTl — 1. 


Then the density of the Ginibre product process with polynomial ensemble initial condi¬ 
tions is given by formula liS. 1\) . 


Proof. We use notation introduced in Section 12.11 To find the density of the Ginibre 
product process with polynomial ensemble initial conditions we need to hnd the density 
of the vector , x™), where x^ = (x},..., xfi are the squared singular values of Gi, 

and x^ = (xfi ..., x^), r G {2,..., mj, are the squared singular values of ... Gi. Note 
the density of = (xj,..., xj,^) is given by equation fl2.2p . To hnd the density of (x^, x^) 
we apply Proposition 14.11 The result is that this density is proportional to 


A(xj,---,x^) 

A(xi,...,xi) 


det 


X, 


2\^2 


(a') 


I', !22 + 1 


e J 


A (xj,...,x^,) det [fiix])]^.^^ [x])'"^ e 


i,j=l i—1 

Continuing in this way we obtain that the density of (x^,..., x”*) is proportional to 


(4.1) 


A (xr, ..., o det [fi{x])]l^^^ det I [xtr {xT") 

X det I ^^rn-2ym-2-yrr.-,-i 


X det (xj) 


1 \ 121 - 1 .^ 2-1 -x^ 

e 3 e 1 



— 2,™.- 
e 3 


*J = 1 


e 3 


... det {xfi 


2\!22-i23-l P. 

' 6 3 


*J = 1 


*4 = 1 


* 4=1 


We compare this formula with equation fl3.ip . The Proposition is proved. 


□ 
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Proposition 14.21 and Proposition 13. ll sav that the Ginibre product process with polyno¬ 
mial ensemble initial conditions is in fact a determinantal point process on {1,, m} x 
M>o. The correlation kernel of this process is dehned by equation fl3.2p . Propositions 
14.3114.61 give the involved transition functions, and the matrix A explicitly. 


Proposition 4.3. Assume that 1 < r < s < m. Then 


(4.2) 


y) 


9r{,^) 1 ^s—r,0 


^r+ly • • • y^s 



y 


where the functions Qr, 1 < r < m, are defined by 


(4.3) 



r = 1, 

r = 2,..., m — 1, 
r = m. 


Proof, li 1 < r < s < m, the two-point function (j)r,s{Xyy) can be written explicitly as 


4^r,s{,Xy 1 /) * ... * (l)s—l,s{,Xy 1 /) 




*r+l 


^s — 1 y 
ts — 2 *s-l 


, \ ^r+1 / ^ \ i^r+2 

^r+1 \ / ^r+2 


X 


^r+1 


y 

ts — l 


dt 


T+l 


tr+1 


0 0 


X 


, 1 = r < s < m, 

-I 

1 < r < s < m, 
1 = r < s = m, 


X 


x'^^e r = 1, s = m. 

To obtain the formula in the statement of the Proposition use the identit}0 


^m,0 

'^0,m 




Xm 

Xo 


(4.4) 


Xo 


hi 


Xi 


b2 


dts-1 

ts — l 


Xr. 


^m—1 


X2 

2^0 XI 


. dxi dXm-l 


Xi 


^m—1 


and the dehnition of the functions Qr (equation fl4.3p h 
Proposition 4.4. The transition functions 

00,s : {1,..., n} X M>o -t C, s G {1,..., m} , 


□ 


^Identity (14.41) can be proved inductively using 




Q-l, . . . , Qjp 

bl,...,bq 


^ \ _ ^m+l,7T, I 0-1; • ■ • ; Op 


V bo,...,bg 


see Luke 
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of the Ginibre product process with polynomial ensemble initial conditions are given by 

fj{y), s = l, 

> . „ / _ \ 

I j (it, s = 2,... ,m. 


(4.5) (/.o,.(j,|/) - <; 


as{y) 


z/2,. . ., z/* 


Proof. The function is the initial transition function for the Ginibre product process 
with polynomial ensemble initial conditions dehned in the statement of Proposition 14.21 
For s G {2,..., m} we have 

OO 

0o,.(j,2/) = 00,1 = / fj{t)(t>i,s{t,y)dt 


0 


dt, 


where we have used the formula for 01^^ derived in Proposition 14.31 

Proposition 4.5. The transition functions 

0r,m+i : IF>o X {1,..., n} -)■ C, r G {1,..., m} , 

of the Ginibre product process with polynomial initial conditions are given by 

e“^a;''i+*“ip(z /2 + t)... T{i'm + i), r = 1, 

(4.6) (pr,m+iix,i) = i x^^+^~^T{ur+i+i) ...T{u^ + i), r = 2,...,m-l, 


□ 


X 


i—1 


r = m. 


Proof, li r = m, then (j)rn,m+i is just the hnal transition function for the Ginibre product 
process dehned in the statement of Proposition 14.21 For r G {1,..., m — 1} we have 


4^r,m+l{x, %) — 0j-,r+l * ... * * 4^m,m+l{x) ^) 


gr{.x) 

X 


OO 



m—r,0 
0,m—r 


0 


^r+1? • • • ? 



OO 

0r,r+l * • • • * (j)jyi—l^jyi(^X ^ t)t dt 

0 


where we have used the formula for 0r,r+i*- • derived in Proposition 14.31 Recall 

that Qr is defined by equation fl4.3p . and that the Mellin transform of a Meijer G-function 
is given by 


OO 



0 


(2i, • • • 5 Ojp 

6l,...,6g 



UT=Xib^ + u)]X=Xa-a^-u) 
ULn+l r(ai + u) ntm+l r(l - - u) ’ 
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see Luke [26|. So 
(4.8) 


^m—r,0 
'^0,m—r 


^r+1 




X 


^dt = a;T(z/j.+i + i)... + i). 


The result follows. 


□ 


Proposition 4.6. For the Ginibre product process with polynomial ensemble initial con¬ 
ditions the matrix A = in formula for the correlation kernel is given 

by 


(4.9) 


hj = T{v 2 +j) ■ ■ ■T{Vm+j) ft{t)e H"'*’ 'dt, 


Proof. The matrix elements aij of A are defined by equation 03.41) in terms of the tran¬ 
sition function 

(f)o,m+i : {1,..., n} X {1,..., n} -)■ C. 

Let us compute this transition function explicitly. We have 

OO 

0O,m-|-l(Lj ) 00,1 *01,m+l(Lj ) ^ 00,1 (b ^)01,m+l (t, J 

0 

OO 

= T{u 2+ j) ■ ■ ■ r(i^m + j) [ fiit)e~H''^^^~^dt. 


Thus Qij is given by equation 04.9p . 


□ 


Applying Proposition 13. II and Propositions 14.2114.61 we see that the correlation kernel of 
the Ginibre product process with polynomial ensemble initial conditions can be written 
as 

Kn,m{r, X] s, y) = x; s, y), 

9s{y) 


where Kn,m{'i^,x', s,y) is defined in the statement of Theorem 12.11 and where the scale 
ia.ctoTgr{x)/gs{y) emerges from the formulas for the transition functions derived in Propo¬ 
sitions It follows from equation fl2.7p that the kernels Kn,m{rj x; s, y) and x; s, 

give the same dynamical correlation functions pi(si,Xi;... ;si,xi), i.e. these kernels are 
equivalent B Theorem 12.11 is proved. □ 


^Two kernels K{x,y) and K'{x,y) are called equivalent if det [K{xi,Xj])F^^ = det [K'{xi,Xj)]F^.^^, 
for any m = 1,2,..., for example K'{x,y) = {f{x)/f{y))K{x,y). Thus two equivalent kernels define 
the same correlation functions. 
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5. Proof of Theorem I2.2L Proposition 12.31 and Proposition 12.41 


In order to prove Theorem 12.21 we use the formulae for the correlation kernel in The¬ 
orem [Til Taking into account that the functions fi{x) are given by equation fl2.8l) . we 
immediately obtain explicit formulae for the entries aij of the matrix T, and for the 
initial transition functions 0o,s(j)S). Namely, 

Ojj = r(z/i-1-i + j - i)r(z/2j)...r(r'mj), e 

and 


(5.1) 


4'o.sij’ y) = a 


s ,0 

0,s 


Ui + j - 1,U2, . . . , Vs 



where j G {1,... , n}, and s G {1, ..., m}. Next we use the fact that the inverse 
of the Hankel matrix 

, hk = r(A; -\- V 

is given bj@ 

V{v + p + l) i-pU-p)i 
^ (r(i/ -h i)fp\ (^1 + ^)k{vi + i)ik\i\ ’ 

see, for example, Akemann and Strahov [7], Section 5. This gives 

\ _ 

^ T{v2 + i) ...T{vm + i) 

^ r(z/i +P + 1) _ 1 

^ (r(i/i -h l))^p! + l)j-i {i - l)!(j - 1)!’ 


where i,j G {!,...,n}. Proposition 13.11 savs that the correlation kernel can be repre¬ 
sented as a sum of two terms. Proposition 14.31 gives the first term for the correlation 
kernel. Denote the second term of the correlation kernel by s,y). The matrix 

A~^ is given explicitly by equation fl5.2p . so we can write 


(5.3) 


where 

(5.4) 


Kn,m{r,x;s,y) 


(j)r,m+l {X, i) {A (j)o^s{j, y) 

i,j=^ 


n—1 


^ ^ ^r,p{x)Qs,p{y) 

p=0 


r(vi+p + i) 

(r(v, + i)fp!’ 


p 

Pr,p{x) ^ ^ 0 r,m-|-l (^5 ^ T 1) 
j=0 


_ {.-v)i _ 

r(z/2 -|- i -|- 1) . . . r(r'm + * + 1)(^1 + l)d! 


^Recall that the Pochhammer symbol, {a)i, is defined by 

{a)i = a{a -I- I)... (a -I- ^ — I). 
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and where 
(5.5) 


Qs,p{y) — '^2 2).j\■ 


The transition fnnction 0r,m+i can be fonnd as in Proposition 14.51 and the transition 
function 0o,s is given by eqnation fl5.ll) . We hnd 

p 

(5.6) 


Pr,p{x) = r(z/i + l)p! ^ 


■1) 


P — 1 


X 


{p - i)!d r(z/i + i + 1)... Ti^Ur + i + 1) ’ 


and 

(5.7) Qs,piy) = T{iyi + l)p\J2-f - -S,-}} I 1 I T 

“ (P-j)!j!r(i/i + 1 + 3 ) 

The residne calcnlations give 


^s,0 

^0,s 


(6.8) 


Pr,p (^) 


T{vi + l)p\ 
2ni 


T{t-p) 


riLo +^3 +1) 


Z/i + j, Z/2 . . . , Vs 


x^dt, 


y 


where Sp is dehned in the same way as in the statement of Theorem 12.21 To hnd a 
contonr integral representation for Qg^p we nse the formnla 


(5.9) G* 


s,0 


+ j, 


fC+ioo 


y]=— T{u + vi+j)T{u + V 2 ) .. .T{u + Vs)y ^du, 

ZilTl —VrvD 


where c > 0. Inserting the contonr integral representation jnst written above into eqna¬ 
tion (15.71) we hnd 


QiM = r(!'i + i)p! 


(5.10) 


^C+200 


X 


2Tii 


r(n + z/2)...r(n + z/,) 


{-1)pPT{u + v^+3) 


+ 1 + j) 


y-'^du. 


The snm inside the integral in eqnation fl5.10p can be rewritten as 

(_l)p-ir(M + + j) _ 1 Y{u + Vi) -A {-p)j{u + Vi)j 1 


(n — V 9 


(p-j)!j!r(i^i + i + j) p! r(i + i^i) (i + «^i)i j! 


E 


i=o 


i=o 


1 r(M vi) 


2+1 {-p,u + Vi; 1 Vi; 1) 


p! r(l vi) 

1 r(MZ/i) (1 - m)p _ 1 T{u + Vi) r(M) 


p! r(l -h vi) (1 -h vi)p p! r(l + vi + p)T{u - p)' 
where we have nsed the formnla (the Chn-Vandermonde snm) 

(c - b)n 


2+1 (-n, b] c; 1) = 


(c)r 
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see Ismail [H], Section 1.4. So 


(5.11) 


Qs,pi.y) 


+1) n;=or(^ + ^,) 

27rir(l + z/i+p) T{u-p) ^ 


Equations fl5.3p . fIS.Sp . and fIS.lip imply that the correlation kernel can be written as in 
the statement of Theorem 12.21 Theorem 12.21 is proved. 

Once Theorem 12.21 is established, Proposition 12.31 and Proposition 12.41 can be proved 
by repetitions of arguments from Kuijlaars and Zhang [23], Section 5. □ 


6. Proof of Theorem 12.61 Proposition 12.71 and Proposition 12.81 

By Theorem 12 .1 1 and Proposition [2]5]the Ginibre product process with initial conditions 
dehned by a truncation of a random unitary matrix is a determinantal point process. 
The correlation kernel of this determinantal process is given implicitly by formula fl2.4p . 
To write the correlation kernel explicitly we need to hnd explicit expressions for the 
matrix A and its inverse. Using equations fl2.5|) and fl2.14p we hnd that in the case under 
considerations the matrix A = is dehned by 

1 

(6.1) aij = P(z/2 + j).. .P(z/m + j) i, j e {1,... ,n} . 

0 

The integral in the equation just written above can be computed using the formula 


1 

0 


P(a + l)P(/j + l) 
P(q( + /? + 2) 


The result is 

ai,j = P(i/2 + j) • • • r(z/^ + 

We need to hnd the inverse of the matrix A = 

, ^ T{i + j-l + ui) 

T{i+j + l-2n)' 

Proposition 6.1. For N = 1,2,and —a, - 


— 1 + i^i)P(/ — 2n — z/i + 1) 
P (2 + j + I — 2n) 

dehned by 
i,j e {l,...,n}. 

-/5 G C \ N the matrix 


f (o + l)i+j \ 

\(q! + /3 + 2)i+jJ 


i,j=0 
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is invertible and its inverse matrix g^ven by 

(_l)*+i (^Q, ^ X). ^ x)^. 


( 6 . 2 ) 

X 




[a + 1)^ [a + 1)^. (a + /3 + 1) 
(2p + a + /3 + 1) (a + l)pp! 


(a + /? + l)p{/3 + l)p{p - i)\i\{p - jy.jl 


——7 {a + (3 + i + 1) (a + (3 + j + 1) . 


Proof. See Theorem 10 in Zhang and Chen 
We apply Proposition 16.11 and hnd 
(6.3) A-i = (7,,):;,.=,, 

where 

1 


□ 


7*,i ~ 


(6.4) 


X 


r(i /2 + *)... r(r'm + i)^{l -2n — Ui + 1) 

T{l — 2n + 2)^ (/— 2?t, + 1) — 2?t, + 1)^._^ 


(- 1 )^ 


T{yi + 1) (z/i + l)j_i(z/i + l)j-i (/ — 2n + 1) 

(2p + I — 2n + 1) {vi + 1) 


pi (I -2n + i) (I -2n + j) 


^ (l-2n + l)p (/ - 2n - z/i + 1)^ (p - i + l)!(i - l)!(p - j + l)!(j - 1)! ‘ 

To compute the transition function 0o,<j(j, y) consider hrst the case when s G {2,..., m}. 
In this case we have 

1 


(6.6) 




We use the formula 

G 


m,n 


(2l 5 . . . 5 Qjp 

6l, . . . 


X 


X 


— 1 \ /^n,m I 1 ^1: • • • ; 1 

see, for example, Luke |25], to rewrite the transition function 0o,s(j:2/) 

I 


( 6 . 6 ) Mj,y) = ‘'"’ 7 ’^ 0 


dt. 


This integral above can be computed, and the result is 

/ — 2n + j 

1/1 + j - 1, z/2,..., Us 

where s G {2,... ,m}, and j G {1, • • • ,n}. To see that the formula for 4>o,s{j,y) just 
written above holds true in the case s = 1 as well use the identity 


(6.7) 


j,»(j. y) = r(; - 2n -1/, + i)g;;' 


s,0 

s 


y 


T{1 -2n-ui + l)Gi;)’ 


1,0 / I -2n + j 


^i+j 
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Now we are ready to derive the formula for the correlation kernel stated in Theorem 12.61 
We start from equation fl2.4p . The second term in equation fl2.4p can be written as 


Kn,m {r, x; s, y) = 


r(Z — 2n + 1) 


( 6 . 8 ) 


V{i>i + l)r(/ — 2n — z/i + 1) 

n—l 


X 


(1 - 


p!(z/i + l)p(2p + / - 2n + 1) 


(/ — ‘In + l)p(^ — 277- — 1^1 + l)p 


(x) QsjP (l/)) 


where 


(-l)P-* {l-2n + l)i{l-2n + i + l)p 


(6.9) P^,p(x) = ^--— , - ii^ Ft 

^ (p - 7)17! (z/i + l)ir(F 2 + 7 + 1) . . . r(z/m + 7 + 1) 


0r,m+l(-^) ^ T 1)) 


and 

( 6 . 10 ) 


Qs,p{y) — ^ ^ 


(l)P-i (/ - 277 + !),■(/ - 277 + J + 1), 


J^O (p-jV-J- 


(z/l + l)j 


iAj + 


Let us derive a contour integral representation for the functions PrAA- have 
0r,m+l(*^7 7 + 1) r(77^_|_]^ + 7 + 1) . . . r(z7^ + 7 + l)x , 

SO equation fl6.9p gives 


( 6 . 11 ) 


r(z/i + i) ^ (-i)p-^ r(/- 2 r 7 + i + i + p) . 

r(z — 277 + 1) (p — 7)17! r(z/i + 7 +1)... T{i>r + 7 +1)^ ’ 


It can be checked by residue calculations that P^ p(x) can be written as a contour integral 


( 6 . 12 ) 


Pr,p{A) 


r(z/i +1) /■ r(f-p)r(t + /- 277 +1+p) ^ 

T{l-2n + l)J lY Tit+ u, + l) 

Sp 


where Sp is a closed contour that encircles 0,1,... ,p once in the positive direction, and 
such that Ret > —1 for t G Sp. 

Next, let us obtain an integral representation for the functions Qs,p(p). The initial 
transition function 0o,s(j + ^,y) is given by equation fl6.7l) . Inserting the expression for 
4>o,sij + l,p) into equation fl6.10l) we hnd 


(6.13) 


- r(z — 277 — z/l + i)r(z/i + 1 ) 

“ r(/- 277 + 1 ) 

(~l)^ ^ r(/ — 277 + j + 1 + p) g 0 / Z — 277 + j + l \ 

^ (p- j)\j\ r(z/i +1 + j) V + j ,V ' 
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The Meijer G-function in the equation above has the following contour integral repre¬ 
sentation 


s,0 


G'l: 


y 


(6.14) 


I — 2n + j + 1 

1 f r(M -F z/i -F j)r(M -F U 2 ) ■ ■ .r(M -h us) 

27ri J T{u + I — 2n + j + 1) ^ 


We insert expression fl6.14p into formula fl6.13p . and interchange the sum and the integral. 
The resulting sum inside the integral can be written as 


(_l)p 3 Y{u + 131 + j)T{l - 2 n+l+p + j) 


(n — i 1i 


^0 ~ -2n + l + u + j) 

= -TTfTT— 1 ^ 1 ] 3^2 (-P, M + 1^1, / - 2u + 1 + p; 1 + z/i, / - 2n + 1 + u; 1) 

T{1 + 13 — l)r(/ — 2n -I- 1 -I- m) p! 

The Pfaff-Saalschiiltz Theorem says that 

3 F 2 (-p, a, b; c, d] 1) = 


(c - a)p{c - b)p 
(c)p(c - a-b)p 


if the balanced condition, c + d = 1 — p + a + b, is satished, see, for example, Ismail [T8] . 
Section 1.4. In our case 

a = u + vi, b = l — 2n + l+p, c = 1 + vi, d = l — 2n + l + u, 
and the balanced condition is satished. Thus we have 

(1 - u)p{i3i -l + 2n-p)p 


3F2 (-p, u + I3i,l - 2n + 1 + p-,1 + I3i,l - 2n + 1 + u-,1) = 


(1 i 3 i)p{ 2 n -u-l-p)p 


Taking into account that 


(1 -u)j 


(w - P\ 


{ 2 n — I — u — p)p {u — 2 n + I + l)p^ 


we obtain the formula 


Qs,p{y) = (- 1 ) 


pr(Z — 2 n — i3i + 1 )T {1 — 2 n + 1 + p)(z/i — I + 2 n — p)p 


r(/ — 2n -|- 1) 


p!(l + vi)^ 


(6.15) 


C+200 


X 


1 f T{u + i3i).. .T{u + i3s) {u - p)^ 


2'Ki 


r(/ — 2n -|- 1 -|- u) {u — 2 n + I + 1 )^ 


-y '^du. 


To write explicitly the second term in equation (12.4p we insert the formulae for Pr,p 
(equation fl6.9p ). and (equation fl6.15p i to equation fl6.8p . This gives the formula for 
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the correlation kernel in the statement of Theorem 12.61 Note also that the second term 
in the formula for the correlation kernel can be written as 


(6.16) 


Kn,m {r,x;s,y) 


1 fUU^iu + iy, + l) 

(27r*)2X__ ^fYyr{t + u, + i) 

S 


n—1 

^(/-2n + 2p + l) 

p=0 


r(t - p)r(t +1 - 2n + p +1) 

r(M + / — 2n + p + 1) 


The sum inside the integral is the same as that in Kuijlaars, Stivigny |23] (see the proof 
of Proposition 4.4 in Kuijlaars, Stivigny |23]), and the rest of the proof of Proposition 12.71 
and Proposition 12.81 is the same as that of Proposition 4.4 and Theorem 4.7 in Kuijlaars, 
Stivigny |23]. □ 


7. Proof of Theorem 12.101 and Proposition 12.111 


When the initial matrix is a Ginibre matrix with a source, the Ginibre product process 
is a determinantal point process on {!,... ,m} x M>o. This follows immediately form 
Theorem 12.11 and Proposition [531 Theorem 12.11 gives a general formula for the correlation 
kernel x', s, y). The hrst term in formula (12.41) for x; s, y) does not depend 

on initial conditions, and it is the same is in the previous cases. Let us derive a formula 
for the second term of Kn^rn{r,x] s,y). By Theorem 12.11 and Proposition 12.91 the second 
term of Kn,m{r, x; s, y) is given by 


Kn,m (r, X] s, y) = ^ ((>r,m+l(T, i) {A . (^o,.(j, 2/), 


*,i=i 


where 

and 

(7.1) (j)o^,{j,y) = 


’,m+i(T, i) P(z/j._|_i T ^)... r(r'm T ) 


oTi (z/i + 1; Qjt) 


r(ni + 1) 


z/2 ..., 


dt. 


The matrix A = is given by 

OO 

ai,j = r(z ^2 + j) • • • r(^'m + j) j oFi (z/i + 1; g^t) 

0 

where i, j G {1,..., n}. Since 


Ll^iy) = 


p^i+ke X; -yt) dt, 


fc!P(z/i + 1) 


0 
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we obtain 




nr(<^o + j)iTi(-9.)e*. 

i=o 

Proposition 7.1. A~^ = This matrix is characterized by 




j=l 1=0 


n 

1=1 

l^k 


-u - qi 

Qk - qi 




Proof. See Forrester and Liu [m Sec. 2 . 1 ], 

Proposition 7.2. Set 

n 

(7.2) ^Pr,m+l{x,j) = ^(t)r,yn+l{x,i){A~^).., 

i=l 

The function (pr,m+i{x, j) can be written as 

OO 

« n 

(7.3) (Pr,m+l{x,j) = / u’'^e~^^r{u]x)e~'^^Y \_~~— 


□ 


t=l 


where 


(7.4) 


"^riu; x) = 


{2'KiyV{i>i + 1 ) 


dwi 


dwr 




71 


7r 


1=1 


X e“i -“’'oF’i z/i + 1; 


XU 


Wi . . .Wr ^ 

In the formula above 71 , ..., 7 ^ are paths starting and ending at —00 and encircling the 
origin once in the positive direction. 


Proof. See Forrester and Liu [T71 Sec. 2 . 1 ]. 

Proposition 7.3. The initial transition function (j)o^s{j,y) can be written as 

yoApy) = ^si-qpy), 

where the auxiliary function (?;?/) is defined by 

C+2CXD 

^s(g;2/) = ^ / + + z). 

c —100 

In the formula above c > —min (z/i,..., Us), s G {1,..., m}, and 

00 

1 


□ 


i9(g;x) = 


r(ni + 1) 


dtpi+z ig 


0 
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Proof. Use equation (EDI. and the integral representation for the Meijer G-function 
inside the integral in equation fl7.ip . 


(7.5) 


1,0 



c+ioo 

f- I + + .-)(!) ‘dz. 

c—200 


□ 


Now we have 


(7.6) 


Kn,m ip 1 T, S, 1/) ^ ^ Pr,m+1 {x, j)0O,s(j) V) 

i=i 




'Ipriu] x) 


n 

E- 

i=i 






1=1 

¥j 


We apply the same argument as in Forrester and Liu [T71 Sec. 2.1], and by the residue 
theorem we obtain a kernel equal to that dehned by equation fl2.18l) . Theorem 12.101 is 
proved. The scaling limits of the correlation kernel in different asymptotic regimes given 
by Proposition 12.91 are derived by a repetition of arguments from Forrester and Liu [13 
Sec. 3.1]. □ 
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